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Abstract 

We introduce a weighted de Rham operator which acts on arbitrary tensor fields 
by considering their structure as r-fold forms. We can thereby define associated 
superpotentials for all tensor fields in all dimensions and, from any of these super- 
potentials, we deduce in a straightforward and natural manner the existence of 2r 
potentials for any tensor field, where r is its form-structure number. By specialising 
this result to symmetric double forms, we are able to obtain a pair of potentials for 
the Riemann tensor, and a single (2, 3)-form potential for the Weyl tensor due to its 
tracelessness. This latter potential is the n-dimensional version of the double dual 
of the classical four dimensional (2, l)-form Lanczos potential. We also introduce a 
new concept of harmonic tensor fields, demonstrate that the new weighted de Rham 
operator has many other desirable properties and, in particular, it is the natural 
operator to use in the Laplace-like equation for the Riemann tensor. 

Mathematics Subject Classifications: 58xxx, 53xxx 
PACS numbers: 02.40.Ky, 02.40.Vh, 04.20.Cv 

Key words: de Rham Laplacian, tensor-valued differential forms, local potentials, 
curvature tensors. 



1 Introduction 

The classical Helmholtz theorem [HE] characterises any three dimensional vector field V 
on Euclidean space in terms of two potentials (scalar and vector respectively), and the 
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proof depends on the existence of a solution to the vector version of Poisson's equation 

o 

V 2 V= V 

o 

for a superpotential vector field V. The potentials follow immediately from the vector 
operator identity 

o o o 

V = V 2 V = V(V- V)-V x (Vx V). 

Four dimensional jSS] and n-dimensional generalisations have been discussed for 

vector fields using analogous arguments. 

More generally, on an n-dimensional pseudo-Riemannian manifold, the (local) Hodge 
decomposition (3JJ (also (23111111312312011201) characterises any p-form £ in terms of a 
pair of potentials (respectively a (p — l)-form and a (p + l)-form), and the proof depends 
on the existence of a solution to the Laplace-like equation 

A £= £ 

o 

for a superpotential S, where A = dS + 5d is the de Rham operator for p-forms [TBI . 
d is the exterior differential and S is the co differential operator. The potentials follow 
immediately (see Result 12.2(1 from this definition 

£ = AS =d(S E) + 5(d S) . 

In the case of general tensor fields on an n-dimensional pseudo-Riemannian manifold 
with Levi-Civita connection V, one is still able to exploit the natural generalisations of 
these three operators by considering the tensor as an appropriate tensor-valued p-form, e.g. 
0I111E!- It is important to note that when a tensor (with more than p antisymmetrical 
indices) is considered as a tensor- valued p-form, the generalised de Rham operator acting 
on such a tensor does not commute with properties involving indices from both its p-form 
part and its tensor- valued part; hence simple natural constructions for a superpotential 
and potentials for a tensor-valued p-form cannot be deduced in an analogous manner as 
for single p-forms. 

Recognising this deficiency, Lichnerowicz ^TJ 02] proposed yet another generalisation 
Al which, when acting on arbitrary tensors, commutes with all their index properties; 
unfortunately, the Lichnerowicz operator Al does not have direct links to the exterior 
differential d and co differential operator 5, and hence to natural potentials in the manner 
of A = d5 + Sd for single p-forms. 

In this paper we overcome all these deficiencies in a different manner by exploiting 
the fact that any tensor can be considered, in a precise and unequivocal way, as an r-fold 
form |49j . and defining a new generalised Laplacian which is in effect a weighted de Rham 
operator A; this operator is directly suited to the r-fold form structure of the given tensor 
and commutes with all index properties, but crucially also has direct links with appropriate 
exterior differential du\ and co differential 8u\ operators (see Section |HJ). Using the new 
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operator A we are able to identify associated superpotentials for any tensor directly via 
a simple Laplace-like equation. Furthermore, we obtain, in an analogous manner, the 
natural generalisation of the Helmholtz and Hodge decomposition: in arbitrary curved 
spaces, when an arbitrary tensor is considered as an r-fold form then it can be expressed 
in a very natural and useful manner in terms of 2r local potentials. By specialising these 
results to double forms, we can then deduce the potential structure of Riemann and Weyl 
curvature tensors respectively, and hence re-establish our recent result which obtained 
one double (2, 3)-form potential for the Weyl tensor in all dimensions [2*7] . 

We begin in Section 2 by formulating familiar results for p-forms in a manner from 
which we can generalise to more general tensors, and in Section 3 we first review how to 
consider any tensor T as an r-fold form Associated with each block, % — 1, 2, . . . r of 
an r-fold form T we have the three generalised operators d^), 5^), A(j) from which we are 
able to determine a superpotential and a pair of potentials; but these superpotentials and 
potentials are highly non-unique with limited practical use. With the introduction of the 
weighted de Rham operator 

A^(A( 1) +A (2) + ...A r ) = i£A w 
' ' i=i 

in Section 4, we establish our basic result that for an r-fold form T there exists an 

o 

associated local superpotential T given by 

A f= T, 

o 

where it is important to note that this superpotential T does not just have the same form 
structure, but also has the same index symmetry and trace properties as T. From this 
we can deduce, in a natural and straightforward manner, a generalisation of the Hodge 
decomposition, demonstrating the existence of 2r local potentials, 

T = AT = l± A (i) T = I E[* W (d W T) + T)\ ■ 

' i=i ' i=i 

This leads naturally to a new definition for harmonic tensors, in general. 

In Section 5 we specialise these results to double (g,p)-forms and identify their four 
different potentials. Further specialisations are made in Section 6: first to (anti) symmetric 
double (p, p)-forms, which are shown to have only a pair of potentials — a double (p, p+ In- 
form and a double (p, p— l)-form — and then to the particular value p = 2. These results 
are immediately applicable to the Riemann curvature tensor — Theorem 16.31 — and to the 
traceless Weyl tensor — Theorem 16.41 [2*7] . 

In Section 7 we demonstrate that the new weighted de Rham operator is the natural 
operator to use in the Laplace-like equation for the Riemann tensor, and we illustrate its 
advantages over the particular de Rham operator which is usually used. In Section 8 we 
summarise, demonstrate how this work provides insights into some other investigations, 
and outline future plans. 
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2 Standard results for £>-forms. 



Let V n be any differentiable n-dimensional pseudo-Riemannian manifold endowed with a 
metric g a b of arbitrary signature. The covariant derivative associated with g is denoted 
either by V or by a semicolon, e.g., V a Vb = v^ a ] and its Riemann and Ricci tensors by 
Robed and R a b = R c ac b, respectively. Our convention for the Riemann tensor follows from 
the Ricci identity: 

2Wa;[fec] = R d abcVd- 

As usual, we use round and square brackets to indicate symmetrization and antisym- 
metrization of indices, respectively. 

The graded algebra of exterior forms is A, with A p denoting the set of exterior p-forms; 
in particular, for p = 0, the elements of A will be scalar functions, (also called sometimes 
0-forms). The canonical volume element n-form r] G A n is denoted by 77 ai ... 0n = T)\ ai ... an ]. 
Then, we can define the standard Hodge dual operator * : A p — > A n ~ p by means of 

SWw^ -^...a^-"* VE6A". (1) 

pi 

It is easy to prove that ** = e(—l) p<Jl ~ p > where e = ±1 =sign(det(g a fc)) is a sign depending 
on the signature. 

One can define a scalar product < , > on each A p in the standard way by first defining 
the function 

(S,$) = S ai ... ap $ ai -^ 
for every E, $ 6 A p and then integrating this over the manifold V n 

< £,$ >= I (£,$)77 = p\ f EA 4 . (2) 

When V n is compact without boundary the integration can be carried out without further 
restrictions; otherwise, at least one of S, $ has to be assumed to have compact support. 
This will be assumed without explicit mention in what follows. This scalar product is 
bi-linear, symmetric and non-degenerate. In the case of proper Riemannian manifolds, it 
is also positive definite. 

The standard operations on the exterior algebra are the exterior differential (also 
called 'curl') d : A p — > A p+l and codifferential (also called 'divergence') 5 : A p — > A p_1 
dH US 123 EH EH • The second is simply defined by 

5 = (-l) p ^d* = e(-l) {n - p){p ~ 1)+1 *d * (3) 

when acting on p-forms. With index notation they can be given as 

(GE) ai ...a p+1 = (P+ l)V[ 0l Sa 2 ...a p+1 ] = (-l) P (p + l)S[a 1 ...a p ;a J1+x ] , 

= —T7 ai y - (—'{\py ; a i 

ai...a v — v < L -'aia2...a p \ ^ ) L - i a^...a v a\ 

for all S G A p . Observe that da = 0, Vet G A n and, as is implicit in this formula and 
follows trivially from the definition (J3J) of 5, for functions / one has: 5f = 0, V/ G A . 
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Elementary properties of these operators are d 2 = and d~ 2 = 0. Moreover, for every 
£ G A p and T G A p+1 we have 

< d£, T > = < E, 5T > 

as can be easily checked by using the Gauss theorem. Therefore, d and S are mutually 
adjoint with respect to < , > jHlE]. As usual, this property can be used to define d and 
5 acting on p-form distributions by means of the standard use of "test p-forms" (i.e., C°° 
p-forms of compact support), see jH]. 

A p-form E is called closed if dH = and exact if E = for some ^ G A p_1 . 
Analogously E is said to be co-closed if <5E = and co-exact if £ = 5T for some T G A p+1 . 
In general, every statement on p-forms has a dual statement replacing d for 5 and the 
form by its Hodge dual. Obviously, all exact p-forms are closed, but the converse is not 
true, globally, in general. The Poincare lemma ensures that every closed form is locally 
exact. Analogously, all co-exact forms are co-closed, and these are locally co-exact. 

The de Rham cohomology class of order p is defined as the quotient of the set of closed 
p-forms by the set of exact p-forms. And the de Rham Laplacian operator A : A p — > A p 
is defined intrinsically by A = dS + 3d (sometimes also written as A = (d + 5) 2 ) jTHj. Its 
formula with index notation is 

(AS) ai ... ap = -V c V c E ai ... ap +p-R c [ ai S c a2 ... ap ] - — ^i? cd [ aia2 E cd a3 ... ap ] . (4) 

Observe that the Riemann tensor terms are not present for 1-forms, while for 0-forms A 
is (minus) the usual Laplacian V C V C = <7 6c V(>V c . 

The following important properties of A are easily obtained (see e.g. [Td1I41[ I2§l l45|) 

Result 2.1 The operator A is linear, self-adjoint with respect to < , > 

< AS,$ > = < E, A$ > VE,$gA p (5) 
and also commutes with *, d and S 

* A = A*, dA = Ad, 5 A = A5 . (6) 
Furthermore, we can immediately prove the identity 

< E, AS > = < dTj, dTj > + < 5E, 5E > . (7) 

In the case of proper Riemannian manifolds this implies that A is a positive operator, 
i.e., < E, AE > > for all E G A p . 

In what follows, and in order to keep the technical complications to a minimum, we 
shall assume that the metric g is analytical. This will allow us to use the simpler versions 
of the existence results for solutions of Laplace-like equations. More powerful results could 
of course be used by assuming C 2 differentiability of the metric, but this will not change 
the main points we wish to make in this paper. Keeping this in mind, a key result for 
p-forms |29[ I4"5| I3U] , which we will generalise in this paper, is 
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Result 2.2 (Local Helmholtz-Hodge decomposition) Given any p-form £ G A p , 

o o 

there always exists a local superpotential S 6 A p such that A £ = £; furthermore there 
always exists a pair of local potentials (\I/, T) with \I/ G A p_1 and T G A p+1 such that 

E = d^ + 5T (8) 

o o 

where ^ — 8 and V = d £. 

Proof. From the structure of the de Rham operator A in (@J), according to the Cauchy- 
Kovalewski theorem [To], we can always find a local solution £, which is another p-form, 

o 

to the equation £ = A E. The remainder of the result then follows from the definition 
A = dd + 6d. □ 

Note that \l/ is missing for the case p = so that in this case the statement is just that 
any function is locally the divergence of a vector field, while for p = n, T is identically 
zero stating that any n-form is closed ergo locally exact. 

An index version of Result 12.21 is 



c 

,;cl 



(AE) ai ... ap = (-l) P p£c[oi...ap-i 5C ap] - (P + l)E [oi ... ap ; £ 
— ( — 1) P P^ r [ai...Oj,_i;o p ] — ^ca\...ap 

= (d^) ai ...a p + (ST) ai ... ap 

= S a ,.. ap (9) 

where we have defined 

o 

^ai.-.ap-i = — E COl .,, 0p _ 1 ' , (10) 
r ai...a p+1 = (-l) P (p+ 1) S[a 1 ...a p ;a p+1 ] • (H) 

Observe that this result is independent of any field equations for E, and is a local result. 
Actually, Result [2~2l can be strengthened in the case of compact without boundary proper 
Riemannian manifolds and one can obtain the global Hodge decomposition theorem pJ71ITol 
El 1111123 113 HHEO]: any £ G A p admits a unique global decomposition as E = d^+SY+T 
where ^> G A p_1 , V G A p+1 and T G A p is a harmonic p-form (that is, T is closed and 
co-closed, i.e., AT = 0) in the same co-homology class as E. Note that the global Hodge 
result can be used to obtain the more modest local Result E21 because, since T is closed 
and co-closed, it is also locally exact and co-exact (see [M]). 

Closely linked to the existence of potentials is their gauge freedom, and there are well 
known results for p-forms which exploit this gauge freedom. However, when one considers 
the more complicated tensor-valued forms which are the subject of this paper, there are 
aspects of the role of gauge which are significantly different from the p-form results. A 
detailed discussion of the role and application of gauge for tensor valued forms will be 
presented in 
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3 Generalization to arbitrary tensors: r-fold forms 



In the previous section the two crucial results were being able to identify a superpotential 
via a Laplace-like equation for the de Rham operator, and then being able to link poten- 
tials in a simple natural manner to derivatives of the superpotential, via this de Rham 
operator, which is of course defined only for differential forms. The question now is how 
to generalize all this to arbitrary tensors. 

Long ago Lichnerowicz |1T) W2\ proposed a generalised Laplacian for arbitrary tensor 
fields which had the first of these properties, as well as a number of other useful properties. 
Given an arbitrary rank-m tensor field T ai ... 0m , the Lichnerowicz operator Al can be 
defined as 

m 

(^L^)ai...o m = — V V c T ai am + R a a T ai ...a s - 1 ca s+ i...a m iX^) 

8=1 

m 

Epc d rp 
- rl - a s at 1 a\...a 3 -\ca a j r \...at-ida t +\...a m ■ 

This operator has a number of very important properties I42j : 

• Al respects the symmetry properties of T ai ___ am ; that is to say, (AiT) ai ... am has 
exactly the same index symmetries as T 0l „ >0m . 

• Al commutes with traces, i.e., the trace on any two indices of (Ax,T) ai ... am equals 
Al applied to the corresponding trace of T ai ... am ; in particular, if T ai ,„ am is traceless 
in any pair of indices, then so is (A2,T) ai ... am . 

• Al is self-adjoint with respect to the scalar product { , } defined by 

{T,S} = [ T ai ... am S a ^ a - rj (13) 

for arbitrary T,S G T m (V n ). (Of course, the same comments as followed (J2J on 
the compactness of V n or of the support of one of the tensors, are in order here.) 
Therefore {A L T, S} = {T,A L S}. 

• In manifolds with a parallel Ricci tensor (V a Rb c = 0, which includes the important 
cases of Einstein spaces and Ricci-flat manifolds,) the following two properties hold; 
see gUHg: 

— when acting on rank-1 tensors, Al commutes with the covariant derivative; 

— when acting on rank-2 tensors, Al commutes with the divergence operator. 

• When acting on p- forms, Al coincides with the de Rham Laplacian A: 

A L S = AS, V S e A p . 
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Unfortunately, the Lichnerowicz Laplacian Al does not have simple natural links to 
first derivative operators for arbitrary tensors in the same way as the de Rham operator 
has for p-forms. 

We believe that the best way to extend the principles of Section 2 to arbitrary tensors 
is to consider tensors as r-fold forms. This terminology was extensively considered in 
and the underlying simple idea is based on the following remark: Given any rank-m 
tensor T ai "' am , there is a minimum natural number r, r < m and a unique set of r natural 
numbers ni, . . . , n r , with Ya=i n % — m i such that T ai " am is a linear map on A ni x . . . x A" r . 
In other words, there always exists a minimum r such that T G A ni g) . . . <g) A n? ., where 
rpai...am j g ^ e appropriate permuted version of T ai " am which selects the natural order for 
the ni, . . . , n r entries. Tensors seen in this way are called r-fold (ni, . . . , n r )-forms [39] . 
In short, all tensors can be considered, in a precise way, as r-fold forms. 

Definition 3.1 (Form-structure number and block ranks) For any tensor T, the 
uniquely defined number r will be called its form-structure number, and each of the Hi, 
the i-th block rank. 

Some simple examples are: any p-form £ is trivially a single (that is, 1-fold) p-form, 
while VS is a double (l,p)-form, with r = 2 and ni = l,n2 = p; the Riemann tensor 
(r = 2, ni = ri2 = 2) is a double (2,2)-form which is symmetric (the pairs can be 
interchanged); the Ricci tensor is a double symmetric (l,l)-form and, in general, any 
completely symmetric rank-r tensor is an r-fold (1,1,. . . ,l)-form. A 3-tensor A abc with 
the property A abc = —A cba is a double (2,l)-form and the corresponding A is clearly given 
by A a bc = A[ a b] c = A ac b. The standard index version of familiar tensors such as Riemann 
tensors R abcd = R[ a b][cd], Weyl tensors C abc d = C [ab][cc[l , torsion tensors T abc = T [ab]c , or 
Lanczos tensors H abc = H[ ab } c , already have the indices in the appropriate permuted 
version, so that they coincide with their tilded versions. In these cases we shall dispense 
with the ~ label. 

The scalar product (J2J) can be immediately generalized to arbitrary tensors of the same 
type (i.e., with the same form-structure number r and block ranks) by first defining 

(T, S) = f ai ... am S ai - am 

and then (same comments as before on compactness) 

<T,S>= I (T,S) V . (14) 

Jv n 

Observe that this scalar product is adapted to the structure as r-fold forms of the tensor 
fields T and S, and therefore it is different from the product defined in (JT3J), so that we 
have, in general, < T,S >= {T, S} ^ {T,S}. This will be relevant in what follows. 
Note, moreover, that { , } is defined for general tensors T, S, not only for those with the 
same form-structure number r and block ranks. As before, the product (|14|) is bi-linear, 
symmetric and non-degenerate. For proper Riemannian manifolds, it is also positive 
definite. 
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Now consider any tensor field T and let r be its form-structure number. For each 
one of the r antisymmetric blocks, one can follow a procedure similar to that recalled 
above for p-forms. To fix ideas, let us select the i-th block, with block rank rii, to treat 
T as a tensor- valued n r form, so that f ai - ah bl _ bni ah + 1 - ak = f ai - ah [ bl _ bn .f h + 1 ~ ak with a 
set of rii < n completely antisymmetrical indices b\ . . . b ni plus a number of extra indices 
denoted by ai ... ah and ah+i . . . a^. These objects are usually referred to as differential 
tensor forms or tensor-valued differential forms, e.g. JHE1IBJE]- Then, the previous 
definition ((TJ) can be trivially extended so that the dual of T with respect to the i-th block 
bi . . . b ni is denoted as *(i)T and defined by jlHl 

(*, ,T\ ai --- a h, , a h+i---a-k = n , rf:a 1 ...a h b 1 ...b n a h+1 ...a k 

As before, *a)*(i) = e(— l) n ^ n_n ^ when acting on T. 

Next we extend the definitions of d, 8, A from the previous section to differential d^), 
co-differential 5^, and de Rham operator Am acting on the i-th block of the r-fold form 
T as: 1 

(rf( i) ^) ai - a % 1 ...^ +1 a ' l+1 - afc = (-l) ni (ni + l)T ai "^ [fel ... bni ^-^ ;bni+l] , (16) 

%) = (-l) n< d W * W = e (-l) (n ^ )K - 1)+1 * W ^)* (i) (17) 
Am = + . (18) 

(An alternative equivalent definition for 5m can be given via dm) and a contraction (trace) 
operator across the i-th and z'-th blocks (HH]; see comments following the definition of the 
trace operator in Section 5.) 

We emphasise that the covariant derivative acts on all indices, i.e., the extra tensor 
indices as well as the explicit form indices. Although the definition for d was extended to 
a tensor valued form in [TJ] and is well known, as far as we are aware all three operators, 
in the form given above, were first introduced in and are less familiar. In this more 
general context, sometimes the extended d is called 'the absolute exterior differential' or 
'the covariant exterior differential'; some references retain d and 5 |E|, some texts replace 
d (5, A) with d (5, A) 0, or with d v (5 V ,A V ) 0IHE2], or with D [HI El EU El El 113 
El El 5 or with V [30 . We shall use the notation of (fTol IT71 ITH|) since we wish to be able to 
combine different de Rham operators. 

The operator dm (5m) produces another tensor with one more (one less) index in 
general, and this new tensor has in general the same form-structure number r. Of course, 
as for single p-forms, there are some special situations: (i) if rii — n then d^T = 0; (ii) 
for any tensor T with form structure number r, VT has form-structure number r + 1 
and can be considered as a definition of 'd( r+ i)T' ; (iii) if rii — 1 then 5mT has r — 1 
as form-structure number. In this case, in order to compute Am one has to allow the 
operator d(j) in the combination d(i)5(j) to act on the missing block — as if 5 m T were an 
r-fold (m, . . . , Tii-i, 0, ni + i . . . n r )-form — in the same way as in (ii). 

1 Whcn applying the operators *m, dm, 8^), we are denning the resulting tensor as having the i-th 
block in the ordered position, so that *(i)T — *(j)T, d^T — dayT, etc. 
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Taking this into account, we need definitions for dmT, 5mT when £ ^ 1, . . . , r, which 
are simply 

d {e) T = VT, 5 {£) T = 0, ££{l,...,r} (19) 

where in the first case the ~ means that the extra index provided by the covariant 
derivative must be placed in the appropriate place within {1, . . . , r + 1}. 

Bearing all this in mind, we can write down explicit index formulas for 5^, A^ and 
A(j) for the cases i G {1, . . . ,r} jf £ : 

(«5 (i) T) 0l - Oh 6l ..^ < _ 1 Bh+1 - a * = — T ai '" ah d >1 ...b n .- 1 ah+1 '" ak ' c > (20) 
(A (£) T) ai ... am = -V c V c T ai ... am (21) 

(A, N T N \ ai --' a h, , a h+i---a-k _ _V7CV7 rp ai ...a h a h+1 ...a k 
K'-^W-l J 6i— 6 Wi ~~ v v c -t 6i...6 nj 

+n l i? 4fel T ai -^ c fe2 ... fenJ ^ 1 - afc - " i(Tl, 2 ~ - - Rcdibi b 2 T ai "' ahCd b 3 ...b ni ] ah+1 '" ak 

k 

-m J2 Rc ai d[bl f -" - 160 ^ 1 "^^..^] ^ 1 - * . (22) 

8=1 

A trivial calculation using the Gauss theorem leads to 

< d (i) T,U > = <T, 5 {i) U > 

where U and T have the same form-structure number r and Ui(T) + 1 = rii(U). Thus, the 
operators d^ and Su) are adjoint to each other with respect to < , > introduced in (fTlj) . 
Observe the similarities and the differences of A/{\ of (|2*2^) with Al as defined in (fT3|) : 

• For each z = 1, . . . , r, Au\T respects the skew-symmetry on the i-th antisymmetric 
block of T, and the symmetries and trace properties on the extra indices not in that 
block; this implies, in particular, that A^T has the same form-structure number 
and block ranks as T. 

• On the other hand, any mixed trace, or mixed index symmetry, involving indices 
from both the explicit z-th antisymmetric block and the rest of the indices is not 
preserved in general. 

• For each i — 1, . . . , r, Auj is self-adjoint with respect to the scalar product (fH|). so 
that < A(aT, S > = < T, A(aS > for all T, S with the same form-structure number 
and block ranks. Moreover, one can prove the identities 

< T, A W T > = < d {i) T, d {i) T > + < 5 {i) T, 5 {i) T > Vie{l,...,r} (23) 

in the same way as (J7J) was obtained. From here, in particular, one deduces that 
A(j) are positive operators in proper Riemannian manifolds: < T, AmT > > 0. 

• When r = 1, that is to say, when acting on p-forms, A^) coincides with the de 
Rham operator: AmS = AE = A^E for all E G A p . In particular, when acting 
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on rank-1 tensors, Am commutes with the covariant derivative in manifolds with a 
parallel Ricci tensor (V a R a b = 0). There is however no simple commutation relation 
between A(j) and S^y, see fffijj) below. 

As with p-forms, we can put forward the following definition: 

Definition 3.2 (i-harmonic and fully harmonic tensors) A tensor field T with form- 
structure number r is said to be i-harmonic, for i G {1, . . . , r}, if and only if Au\T = 0. 
Such a tensor will be called fully harmonic if it is i-harmonic for all % = 1, . . . , r. 

Note that the harmonic property in the sense of Lichnerowicz (i.e., A^T = 0) is different 
from these new harmonic properties. This may have important consequences for the the- 
ory of harmonic tensors in Riemannian or pseudo-Riemannian manifolds. See Definition 
14. II and the results which follow it. 

It is important to note, unlike the situation for p-forms, that d?U ^ and 5fy ^ in 
curved spaces. Using the Ricci identity, one obtains 

(dli)T) ai '' Mh b 1 ..K i +2 ah+1 ''' ak 

= hm + l)(n, + 2) £ R a ° c[bni+lKi+2 f^^^^j^^ , (24) 

Z s=l 

1 k 

{ X 2 T\ a i--- a h a,h+i—ak pa s de rpa 1 ...a s - 1 ca s+ i...a h a h+1 ...a k (OK\ 

I 6i...6 ni _ 2 — 7^Z^ n c 1 debi...b n .- 2 ■ \ zo ) 

Z s=l 

Note that for higher derivatives all derivatives of the Riemann tensor disappear because of 
the Bianchi identities; this means that even derivatives will give terms involving products 
of the Riemann tensor and T, while odd derivatives will give terms involving products of 
the Riemann tensor and a derivative of T. 

have a similar product structure with the Riemann tensor, since all derivatives of the 
Riemann tensor disappear because of the Bianchi identities. 

Observe that, in flat spaces, these operators are nilpotent. (Compare with [TH], where 
similar ideas are developed in flat Euclidean spaces.) Note also that = whenever 
ni G {n — 1, n}, and that $~ = if rij = 1. More generally, = for m > n — rij, and 
5™\ = for m > Hi, so these curved space operators are also nilpotent, but depending on 
the dimension of the space and the block rank of the tensor-valued form being acted on. 

The commutation properties with A(») follow directly from ()24H25J) using 

d(i)A (i) - A(j)d(j) = d(j)5(j) - %)d(j) , %>A(i) - A (i) 5 (i ) = S^d^ - d^Sf^ . (26) 

(Observe that, in flat spaces, these operators commute.) Finally, A(j) commute with each 
of the Hodge dual operators *q) defined in ffTB]) 

A « = A W Vi, j G {1, . . . ,r} . 

The usefulness of these operators (fTB^). (fTTj). (JTSJ) is that we can now translate the 
calculation performed in Q to arbitrary tensors. Although the operators are different 
from the de Rham operator A it is still possible to apply the Cauchy-Kovalewski theorem 
[To] , and we can generalise Result 12.21 to 
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Result 3.1 Let T be any tensor field, and let r be its form-structure number and (n±, . . . , n r 

o 

its block ranks. There always exist local superpotentials T with the same form- structure 

o 

number and block ranks as T such that A^ T— T; furthermore, for each i = 1, . . . , r 
there always exist a pair of local potentials (Y^,Z^), such that 

o o 

T = T) + <%)(%) T) = 8(i)Yw + d(i)Z(i) 

o o 

where we have defined the r-fold form potentials = d^ T and = 8^ T. 
In index notation: 

o _ 
( A,^T\ ai --- ah , , a h+l---a-k — rpax...a h a h+1 ...a k 
V W / b 1 ...b„ i — J- bi...b ni 

and 

rf>a\...a h a h+l ...a k _ ^a x ...a h a h+1 ...a k ;c, /_i \rn-l ya\...a h a h+l ...a k (07\ 

1 bi...b ni — ~ Y cb!...b ni +\ L ) n i L [bi...6 nj -i ;6nJ 

where Y ai - ah bl .^ bn , +1 ah + 1 ~ ak is a rank-(m + l) tensor of type r-fold (n 1? . . . , + 1, . . . ,n r )- 
form and Z ai - ah bl ^. bn _ 1 ah+1 - ak is a rank-(m — 1) tensor of type r-fold (ni,...,7ij — 
1, . . . , n r )-form, given respectively by 

o 

■yai...a h a h+1 ...a k _ /_i \raj / . , i \rpa\...a h a h+1 ...a k 

1 6i...6„ i+ i — \ 1 ) \'h ~ L ) ± [bi-b ni ;b ni +i] ? 

Z ai - ah bl ... bni _, = -fo 1 "- ah cb 1 ...b n .- 1 ah+1 " Mk ' c (28) 

Trivially, these potentials in the above result are highly non-unique. For r-fold forms 
with r > 1 we can repeat the above constructions for each % e {1, . . . ,r}. Therefore the 
previous definition of potentials can be applied to each of the antisymmetric blocks, and 
thus there are many possible different pairs of potentials. 

A second disadvantage is that each of the operators A^, i = l,2,...r, in general, 
does not have all of the useful properties which we would like. Each operator is adapted 
to respect the index properties of its corresponding block, but not for the other blocks of 
indices. So, for instance, any index symmetry or trace property across the different blocks 
of a tensor will not necessarily be directly reflected in its superpotential. 



4 A weighted de Rham operator and associated po- 
tentials. Harmonic tensors. 

Of course, one can mix the operators A^, i = 1, 2, ... r, weighting them, and obtain new 
Laplace-type operators with similar or better properties. A particularly good one is, 

Theorem 4.1 The operator A given by 

A = ;( A (D + A ( 2 ) + • • • + A (r) ) = I E A (0 (29) 

' ' 1=1 
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is linear, self- adjoint with respect to the scalar product \14\l , respects all index symmetry 
properties, and commutes with all trace operations when acting on arbitrary tensor fields 
(here r is the form- structure number of the tensor field). It is related to the Lichnerowicz 
operator by 

A = ^A L - ^V C V C . (30) 

Proof. The linearity and self-adjointness follows from these properties for each of the 
A(j). The formula 1)30)1 can be easily established from (jl3j) and (J22j) (see also formulas 
(134)1 and (J40M4)) below). The commutativity with index permutations and traces then 
follows immediately from the properties of and V C V C . D 
So we find that although A L and A do not coincide, they are closely related, and 
indeed A has most of the useful properties of A L ; but crucially A in addition has direct 
links with and 5^. Hence we believe that A is a more powerful alternative than the 
classical Lichnerowicz operator A^ since it is so well adapted to dealing with the r-fold 
form structure of the tensors. Of course, for single p-forms we have 

AS = AS = A L S, VSgA p . 

An extremely important consequence of this new operator A defined in (129)1 is that for 

o 

any given tensor field T, there exists an associated superpotential T, by which we mean a 
superpotential, not just with the same form-structure number and block ranks, but also 
with the same index symmetries and trace properties as T. Moreover, from this associated 
superpotential, a set of 2r potentials can be obtained in a natural and straightforward 
manner from the definition of A 

Af = i£[* (i) (d (0 T) T)\- 
r i=i 

Explicitly we have thus proven, 

Theorem 4.2 Given any tensor field T , there always exists an associated superpotential 

o _ o 

T such that A T— T; furthermore, if r is the form- structure number of T, there always 
exists a set o/2r local potentials (Ym, Zu\), i = 1,2, . . .r, such that 

1 r 
' i=l 

o o 

where = T and Z^) = T are the potentials. 

In our earlier derivation of the potential for the Weyl tensor [27j the operator which we 
constructed for the superpotential is the special case of A for a double (2, 2)-form. To the 
best of our knowledge, this seems to be the only time such a weighted de Rham operator 
A has appeared in the literature. There are a variety of other generalised Laplacians 
[3 E21 EH] but for tensor- valued forms in pseudo-Riemannian spaces with Levi-Civita 
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connections the only Laplacians mentioned are usually — V C V C , Lichnerowicz's Laplacian 
Al and the Laplacians which we have labeled A^. 2 In particular for double (2, 2)-forms 
such as the Riemann tensor, the Laplacian used in 0111111113111, is the one we label A( 2 ), 
and as we argue in Section 7, this is not really a suitable operator in that context. It is 
only in ^1] that any attempt is made to compare different Laplacians, where it is pointed 
out that, for the special case of a symmetric 2-tensor, the Lichnerowicz Laplacian can be 
written as A^ = Am + A( 2 \ + V C V C (which can be considered as a special case of formula 
(jnOJ)); but the significance of (A^ + A( 2 )) as a Laplacian was not recognised in [Tlj . 

We start with a simple logical definition which will allow us to show the potentialities 
of the new Laplacian A with a very simple preliminary application. 

Definition 4.1 A tensor field T will be called harmonic if and only if 

AT = 0. 

Obviously, any fully harmonic tensor in the sense of Definition 13.21 is trivially harmonic. 
The converse, however, does not hold in general. Nevertheless, from formula (J23j) we 
deduce the general identity 

< T, AT > = - < T, A ( j)T >=-£(< d (*) T > rf « T > + < %) T > h) T >) ( 31 ) 

r i=l r »=l 

for arbitrary tensor fields. Therefore, it is straightforward to obtain the following converse 
in proper Riemannian manifolds. 

Theorem 4.3 Let (V n ,g) be a compact without boundary proper Riemannian manifold. 
Then, a tensor T is harmonic if and only if it is fully harmonic, and if and only if 

d (i) T = 0, % } T = 0, Vz G {1, . . . , r} . 

Remarks 

• The above definition is different from the Lichnerowicz harmonic property which 
makes use of the condition A^T = 0; and we emphasise that, in general, there is no 
analogue of this theorem by using A^. 

• This definition differs from the harmonic definition in for 2-fold forms in proper 
Riemannian space; in fact the definition there coincides with our definition of 1- 
harmonic 2-fold forms, which is equivalent to Aiu = where a; is a (p, g)-form. 

• See Section [7| for details concerning the harmonic property of the Riemann tensor. 

2 To avoid any misunderstanding we point out that what we have called the 'usual Laplacian' — V C V C 
(and with the symbol Am in is sometimes given by the symbol A, 021 E3 and called the 'rough 

Laplacian' [121 or the 'naive Laplacian' |14j . 
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• The tensor decomposition into superpotentials and potentials for vectors and 2- 
tensors developed in [32] exploited the Lichnerowicz operator A^ in order to find 
the superpotential U for any tensor V via A L U = V; however, such a superpotential 
U does not, in general, yield potentials in a straightforward manner as there is no 
direct natural link to any first order operators. 

For general calculations, we would need the commutation properties of the operators 
d({\ and 5(j). Actually the general case requires a detailed and cumbersome study in non- 
flat backgrounds; but we are only interested in the case of double forms (r = 2) in this 
paper, and so we devote the next section to this case in full detail. 



5 Double p)-forms 

In this section we deal with tensors with form-structure number r = 2 and block ranks 
Hi — Q an d n 2 = p, called double (q, p)-forms. These are tensors with p + q indices 
which fall into two categories: a block of q antisymmetric indices, and another block of 
p antisymmetric indices. By rearranging if necessary, the tilded version of the tensor can 
then always be written as 

rp ai ...a q _ if>[ai..Mq] 

J- bi...bp — -L [&1...&J,]- 

In addition to the Lichnerowicz operator, there are three other Laplace-like operators 
A(i), A( 2 ) and A as defined in the previous sections acting on these objects. Using (fT%|). 
and (|29p their explicit formulas are 



(A ( i)T) ai -%,., p = -V c V c T ai -%,.. bp + qB c[a "T c a *- a \... bp (32) 

Q(q- 1 ) R cd[a 1 a 2 f cd a 3 ...a q ] bi ^ + qp cd[bi f a ^ cd b2 ... bp] ; 



(A( 2 )T) ai -° 9 6l ... 6p = -V c V c T ai - a ' b ,.. bp + pR c[bl f a ^ c b2 ... bp] (33) 

-^% t2 r-n 3 ...y+^ [fll +^'''"' ] ^,..y ; 



(AT) s - a %,., p = _v c v c T ai - a v.. bp + ^R c[ai f c "^ bl ... bp + ^R c [ b J"« c b2 ... bp] 

+qpR [ai c d[b J"^ d b2 ... bp] (34) 

The same comments as before for the cases with p or q equal to 1, n — 1 or n are in order 
here. We also note from formula (j3*Uj) in Theorem l4.1l that A is related to the Lichnerowicz 
operator in this case by 

A = i(A L -V c V c ). 
Specialising Theorem 14.21 we obtain, 
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Corollary 5.1 Given any tensor field T with the structure of a double (q,p)-form there 

o _ o 

always exists an associated local superpotential T such that A T— T; furthermore there 
always exist local potentials Y"(i),Y"( 2 ), Z(i), Z^) such that 



T =2 ( s o-) Y w + 5 m Y m + d a) z (i) + d m z m) > ( 35 ) 

o o 

where the double (q + l,p)-form = <i(i) T, the double {q,p + l)-form Ym) = dp) T, the 

o o 



double (q — l,p)-form Zm = <5(i) T and the double (q,p — l)-form Zm = 6(2) T . 

In the sequel we will need to deal with contractions as well as permutations of indices 
for different expressions, so we make two formal definitions which will be useful. 

Definition 5.1 (Trace of a double form) The trace of a double (q,p)-form T is the 
double (g — l,p — l)-form tr(T) given by 

( t r7r)r-% 2 ... 6p = f ca ^ cb2 ... bp . (36) 

Once more, we must remark that if q — 1 (or p — 1), then the first (second) block 
disappears after taking the trace tr, so that the resulting tensor has a form-structure 
number less than 2. In these situations, and as happened in other situations explained 
before, sometimes it is necessary to consider the resulting tensor tr(T) as an equivalent 
double (0,p — l)-form (or double (q — 1, 0)-form). Bearing this in mind, some useful results 
follow immediately from the respective definitions (|TB|) and (120}) . 

tr(d ( i)T) = -d(i)tr(T) - 5 (2 )T, tr(rf( 2 )T) = -d (2 )tr(T) - <5 ( i)T, (37) 

tr(5 (1) T) = -5 ( i)tr(T) (q > 2); tr(5 (2) T) = -5 (2) tr(T) (p > 2). (38) 

In [23 the second of equations (|3Tjl was used as the definition of for a double (p, q)- 
form. 

Definition 5.2 (Transpose of a double form) For a double (q,p)-form T we define 
the (generalised) transpose *T of T as the double (p,q)-form given by interchange of the 
blocks: 

(t r r\ a i--- a p — nh a 1 ...a p 

\ -L ) b 1 ...b q — -L b 1 ...b q 

Obviously U T = T, and the following useful results also follow from the definitions 

d {2) ( t T)= t (d {1) T), 5 {2) CT)= 4 (5 (1) T) (39) 
In addition, we note that the operations of transpose and trace commute, 

4 (tr(T)) = tr(*T) . 

Of course, we could also define traces and transposes of the more general r-fold forms 
in the previous section, by taking pairs of blocks at a time. 
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For future reference we give explicitly the five cases of (1, 1)-, (2, 1)-, (2, 2)-, (2, 3)- 
and (2, 4)-forms respectively for (JMj) 

- (AT) ai 6l = V c V c T a \ - ^R ca >T cbl - l -R cbl T a ' c - R a \ dbl T cd (40) 

- (AT) aia2 fel = V c V c T aia2 6l - R c[ai f c a2] bl + ^R cdaia2 f cdbl (41) 



(AT) aia V 2 = V c V c f aia2 blb2 - R c[ai f c a2] blb2 + R cda ^f ( 



^ 1 yam \ m:- " / 

cdb\b2 



-R c[bl f a ^ c b2] + X -R cdWb2 f a ^ cd - 4R^ cd[bl f a ^ cd b2 ] (42) 

/ \rp\a\a,2 V7 C V7 r f'aia2 T? c \ a i r F a 2] i ^ T>cda\0i2rp 

— y^ 1 ) 6ifef*3 — V V c i — -K J-c b!b 2 b 3 + J-cdb^bz 

3 ~ 3 
~ 2^c[b 1 T aia2C b2bi \ + -R cd [ blb2 T aia2Cd b , A ] — 6i?' ai cd [b 1 T a ^ cd b2bz ] (43) 

/AT^ ! ^ V7 C V7 r pa\a2 r>c\a\rp 0,2} 1 ^ Tfcdaxa^rn 

— y^ 1 ) &if> 2 f>3&4 _ V V c J b!b2b 3 b 4 — -K i c J 6 1 6 2 fo 3 fo 4 + ~rC J- cdb 1 b 2 b 3 b4 

— 1 ^Rc[b 1 T aia2C b2b3b4 ] + 3-R c d[;, 1 ; )2 T aia2Cd b36 4 ] — 8R^ ai c d [b 1 T a2 ^ cd b2b3b4 ] (44) 

For the case of double (q, p)-forms, a straightforward computation provides the commu- 
tation properties of the operators and 5(j) for i,j G {1,2} 

([rf (1) ,d (2) ]T) ai - a9+ \ 1 ... fep+1 = L_i_(p+i)( g + i) x 

(g^c[ 6p+1 Ka9+i r ai - 09 - llc 6l ... 6p] - ^ cK+1 [ bpf , p+1 T ai -^ 6l ... fcp _ l]c ) , (45) 

(Ki),^)]^" 1 -^ 1 ^..^, = (-l)*(g + l) (|^ [a ' a ' +1 T--^- 1 ]V... fep _ 1 

£(1),£(2)]^) 9 6L..6P-! = ^"2— R Ce d[biT dai - aq - 1 b 2 ...b p - 1 ]ce 

-^Rj^f a2 --- a ^ ce dbl ... bp ^ (47) 

(Observe again, that in flat space these operators commute.) In will be useful in subse- 
quent calculations to note from (|3~3"j) and 

(A(i)T)ai ...a v ^ _ (A (2) T) ai - v.. 6p = qR c[ai f c a2 - a ^ bl ...b p - pRc [bl T ai - a «%... bp] (48) 

?(? ~ 1) D cd[aia 2 ^ a 3 ...a„l , P(P ~ 1) D 
^ -K J cd V..&P H ~ -Kcd[fei6 2 J 6 3 -6p] " 
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6 Double (p,p)-forms: Curvature tensors 



The transpose *T of T is of special relevance for the special case of the double (q, p)-forms 
in which both blocks have the same number of indices, say p — these double (p,p)-forms 
include, of course, the important case of curvature tensors. In this case the two blocks 
of T can be interchanged, and hence we will say that a double (p,p)-form is symmetric 
if T = *T, and antisymmetric if T = — l T (in this case only for p > 1). Of course, 
for p > 1 any double (p,p)-form can be decomposed uniquely into a symmetric and an 
antisymmetric one (if p = 1 the decomposition gives a 2-form and a double symmetric 
(1, l)-form). Hence, without loss of generality, in what follows we will only consider these 
two cases, 

T = ±*T or f ai -°V.. 6p = ±T bl ... bp ai '" a * . 
Then it follows trivially from that 

d m T = ± \d {x) T), 5 {2) T = ± \5 {1) T) (49) 

and therefore A( 2 )T = ±*(A(i)T) so that 

AT=i(A (2) T±'(A (2) T)). 

o 

Since the associated superpotential T will have the same symmetry properties as T, it 
follows that the four potentials Y(i), Y(2), ^(i) and Z(2) defined in the previous section 
satisfy 

y (2) = ±*y (1) = (so) 

Z {2) = i'Zpj^C-l^Z, (51) 

where furthermore the completely antisymmetric part v4[Y( 2 )] of Yt 2 ) (or Y) may vanish 
identically due to the (anti) symmetry of the associated superpotential: 

A[Y {2) ] = for { T = tT tT P ° dd ' (52) 
I 1 — — 1 p even. v ' 

The index notation for this restriction is simply 

. j T = l T p odd, , s 

y [ a 1 ...a p 6 1 ...6 p+1 ] = for | T = ^ ^ (53) 

Theorem 6.1 Given any tensor T with the structure of a double (anti) symmetric (p,p)- 

o _ o 

form there always exists an associated local superpotential T such that A T— T; further- 
more there always exist a pair of local potentials Y( 2 ), Z( 2 ) satisfying A5ty) such that 

T= l - [5 (2) Y {2) ± % 2) Y i2) ) + d (2) Z {i) ± \d {2) Z {2) )} , (54) 
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where the double (p,p + l)-form Y( 2 ) = d(2) T , and the double (p,p — l)-form Z( 2 ) = <5(2) T. 
Using the notation introduced in A5U\) and \51)) . \54\ l can be written in index notation 

as 

T ai - a V.. bp = ^(^ ai - ap 6 ,..6 pC ;C±y 6i-^ ai "' ap % 

+ ^ ai - ap [; > 1 ... 6p _ 1 ; ; > p ]±^ 1 ... fep [ai - aj '- i;ap] ) (55) 
where the potential Y satisfies A5ty) . and the potentials themselves can be given in terms 

o 

of the double (anti) symmetric associated superpotential T by 

o o 
-i^ai...a p _ / , -i \rpai...a p ryax...a p _ rpa x ...a p ;c 

I b!...b p + 1 — —{P + J-JJ P [b 1 ...b p ;b p + 1 ], * bl...bp-l -L Cbl.-.&p-l ■ 

Since the operations of transpose and trace commute, then if T is (anti)symmetric, so 
is tr(T). Suppose now that, in addition to the (anti) symmetry between blocks, the double 
(p,p)-form T is traceless, i.e., 

tr(T) = . (56) 

o 

The associated superpotential T will also be traceless with 

tr(T) = . (57) 

But then the potentials Yp) and Z( 2 ) of (|54|) are not independent and one can easily check 
using (JSTj) that Z^) is essentially the trace of the transpose of Y( 2 y. 

Z {2) = T tr('F (2) ) 

or with indices and the notation of (|50H5H) 

Z ai -^ bl ...b p ^ = TY cbl ... bp _ 1 cai - ap ■ (58) 

Furthermore 

tr(tr(F (2) )) = 

which becomes in index notation 

^6 1 ... fep _ 2 cdai - ap - 1 =0. (59) 

Therefore, we have 

Theorem 6.2 Given any tensor T with the structure of a double (anti) symmetric trace- 

o _ o 

less (p,p)-form there always exists an associated local superpotential T such that A T— T; 
furthermore, there always exists a double (p,p+l)-form local potential Ym satisfying KUty) 
and tr(tr(Y(2))) = such that 

T = \ (5 {2) Y (2) ± % 2) Y (2) ) - d (1) tr(y (a) ) =F *(d (1) tr(y (2) ))) , (60) 
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where the double (p,p + l)-form Y^) = ^(2) T. 

With the notation introduced in A5(J\) . the index version of A6'U\) can be written as 

rpa\...a v _ _ (■yai...a p \ c -\- V a\...OpC 

1 bi...b p 2 \ bi...bpC - n 1 b 1 ...b p ;c 

-pY c[a ^-\ bl ...b p M T P Y c[bl ... bp . 1 cai - ap ;b p] ) (61) 
where the potentialY satisfies i f55j) and KUty) . and is given in terms of the double ( anti) symmetric 

o 

traceless associated superpotential T by 

Y a ^ bl ... bp+1 = -(p + l)f ai -° p [6 1 ...6 pi 6 p+1 ] • (62) 

In j2B] it was shown in 2m dimensions that a traceless symmetric (m, m)-form has a 
traceless (m, m — l)-form potential. That result is a special case of the much more general 
result in this theorem, where the double dual of the (m, m — l)-form potential in is 
precisely the (m, m + l)-form potential Y in this theorem. 

For future reference we give explicitly the index version of the potential structure for 
a symmetric 2-tensor 

T ab = \{Y abc ^ + Y bac ^ + Z a , h + Z b , a ). (63) 
and for a traceless symmetric 2-tensor 

Tab = - iXabc C + Ybac' C ~ Y C ca . h — Y C cb . a ) . (64) 

where (in both cases) the double (1, 2)-form Y a bc satisfies Y[ a 6 C ] = 0. 

There is a well known decomposition for symmetric 2-tensors in three dimensional 
spaces, but it is restricred to proper Riemannian space ^7J E3] , and so differs from this 
one. Of course, the antisymmetric 2-tensor is just a single 2-form whose decomposition 
(JSJ) is easily seen to agree with the decomposition which follows from Theorem 16.11 

6.1 Application to general curvature tensors 

Let us apply the above results to the case of Riemann candidates, that is, tensors with 
the algebraic properties of a Riemann curvature tensor. Let 7Z a bcd be any such Riemann 
candidate, that is to say 

Ttabcd = Tt[ab}[cd\-> T^a[bcd] = (=>■ 7Z abc d = Tlcdab) j (65) 

so that TZ abc d is in particular a symmetric double (2,2)-form. 

o 

Let TZabcd be the associated local superpotential for 1Z ab cd by 

_ o 

{^7t)abcd = TZabcd 



(66) 
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and because of the properties of A, 1Z is also a Riemann candidate; the pair of potentials 
defined in Theorem 16. II become now 

Yabcde = y[a6][cde] = — 37Z ab [ c d ]e ] , Z abc = Z[ab] c = —1Z ab dc ■ (67) 

o 

It is important to realize, since 1Z is a Riemann candidate satisfying (J65|) . that these 
potentials satisfy the additional symmetries 

Y a [bcde] = 0, Z[ abc ] = 0, (68) 

the first of which implies the following useful properties 

y C [bcd]e = 0, yjaficdje = 0, Yabcde = 3Y[cde]ab = 3Ya[cde]b, Y a [bc]de = —Ya[de]bc- (69) 

o 

Given that (|66j) always has local solutions 1Z for any given 1Z we have, 

Theorem 6.3 Any Riemann candidate tensor lZ ab c d has a pair of local potentials given 
by a double (2, 3) -form Y ab cc i e and a double (2,l)-form Z ab c with the properties $6'8\) such 
that 

Ttabcd — — (Y a bcde' e + Y^abe^ + 2Z ab [ c . d ] + 2Z cd [ a . b ]j . (70) 

The potentials themselves can be given in terms of the associated (Riemann candidate) 
local superpotential lZabcd by ffify . 

Of course, any Riemann candidate can be decomposed in terms of its trace (its 'Ricci 
tensor') and its traceless part (its 'Weyl tensor'). The trace is a double symmetric (1,1)- 
form, and therefore the general Theorem 16.11 applies, and in particular ()63j) . After a little 
rearranging, it can be seen that the potentials for this Ricci part are essentially the traces 
of the pair of potentials defined in (J67|) and the formula relating them is the trace of (J7Uj) . 



Corollary 6.1 Any Ricci candidate tensor lZ ab has a pair of local potentials given by a 
double (1,2) -form Y c ab with the property Y[ abc ] = and a double (l,0)-form Z a such that 

TZ ab = X - (Y abe ' e + Y bae ' e + Z a , b + Z b , a ) . (71) 
The relation with the Riemann candidate local potentials of J7fl| ) is 

yc yec y c y y e 

1 ab 1 eab ^ab t ^ea • 

The traceless part of the Ricci candidate tensor (essentially any traceless symmetric 2- 
tensor) has an even simpler structure requiring only one potential, as shown in (|64p: 
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Corollary 6.2 Any traceless Ricci candidate tensor TZ a b = !Z a b — 9ah^, where 1Z = TZ C C , 
has a local potential given by a double (1,2) -form Y c a b with the property Y[ a bc] = such 
that ^ 

TZab = ^ (Yabe' e + Y bae ' e ~ Y C ca - h ~ Y C cb;a) ■ (72) 

The relation with the Ricci candidate local potentials of \71\j is 

2 2 
Y c a b = Y c ab H — g c [ a Y e b] e H — 9c[ a Zb] ■ 
n n 

More interesting is the traceless part of a Riemann candidate, which we call a Weyl 
candidate, that is, a double (2,2)-form C a b c d = C[ Q 6]cd — Cab[cd[ with the algebraic properties 
of the Weyl conformal curvature tensor: 

C a bca = 0, C a [bcd] = 0, (=>■ C a b c d = Ccdab) , (73) 

so that C a ft c d is in particular a traceless and symmetric double (2,2)-form. By considering 

o _ o 

the associated superpotential C for C given by C = AC, we can then deduce from Theorem 
16.21 and Theorem 16. H| 

Theorem 6.4 Any Weyl candidate tensor field C a b c d has a double (2, 3) -form local poten- 
tial Pabcde with the properties 

Pa[bcde] = 0, P^abc = (74) 

such that 

= \ + iV* « - - 2^1^/1) . (75) 

The potential itself can be given in terms of an associated local (Weyl candidate) super- 

o 

potential C ab c d by 

P ab cd e = P [ab] [c d e] = -SC ab [cd]e] • (76) 

This result was originally obtained, by a more direct route, in j2Zj. The operator A 
which we have used in this theorem, specialising via (J42|) . is easily seen to coincide with 
the operator which we constructed for the superpotential of the Weyl candidate tensor 
in |2Zj. As noted there, in four dimensions our new (2, 3) -form potential P ab c de coincides 
with the double dual of the classical Lanczos (2,l)-form potential H ab c for the Weyl tensor 

Immediate consequences from the first of (|74*j) are the following useful properties, 
reminiscent of (|69j) 

P e [bcd]e = 0, P[abcd]e = 0, P^ ode = 3P[ c de] ab = 3P^ a [ c d e ]^ , P a [bc]de = ~ Pa[de]bc ■ (77) 

The basic equation (|75|) can be written in several alternative but equivalent forms by 
using properties (|74|) and (|77j). Other possibilities are 

Cab r>ab ;e , ryelab] r> abe TD^W 'M /' r 70^ 

cd — r c de + f ' c d;e ~ *e[c ;d] ~ " cde {'o) 

and 

C ab cd = P^cde^ + P^cd^ ~ P^cde^ ~ 2P^ e [c;d]. (79) 
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7 Laplace-type equation for the Riemann tensor 



The intrinsic, as well as quickest and simplest, way to obtain a Laplace-like equation (a 
wave equation in Lorentzian signature) for the electromagnetic field F is as follows. The 
2-form F satisfies the Maxwell equations 

dF = 0, 5F = J (80) 

where J is the electric current 1-form. By using the de Rham operator the Laplace-like 
equation follows from these equations strightforwardly 

AF = dJ. (81) 

As is known, in the absence of electric charges and currents (J = 0) one has the simpler 
equation 

J = =>■ AF = 0. 

The index version of (|81|) is the well-known equation ("wave" equation in Lorentzian 
signature) 

V c V c F afc + R efab F ef - R ae F e b + R be F e a = J [a ' b] 

or equivalently 

V c V c F ab - 2R eafb F ef - R ae F b + R be F e a = J [a;6] . 

As has been always the case, the electromagnetic field may serve as source of inspiration 
for the gravitational field, and thereby for the curvature tensor of any pseudo-Riemannian 
manifold. Thus, there have been some attempts to formulate the analogue of (|81|) for the 
Riemann tensor with a generalised de Rham operator, using the Bianchi identities as the 
analogue to the Maxwell equations (|80p. We will use 3? to represent the double symmetric 
(2,2)-form defined by the Riemann tensor R ab c d- Then, the Bianchi identities R a b\cd;e] — 
can be written in a concise form using the notation of this paper as 

<f( 2 )£ = 0. (82) 

Hence, using ()37|) to take the trace of (J%2~)) we readily obtain 

5 {1) 3t = -d {2) mc (= <f( 3 )tr(£)) (83) 

where 9fac = tr(3?) denotes the double symmetric (l,l)-form representing the Ricci tensor 
R a b. (Alternatively ()83|) could be obtained by translating from the contracted index 
version of the Bianchi identities R a bcd' a = —^Rb[c-,d}-) 

The symmetry of the Riemann tensor = 3? implies that (J82)) is actually equivalent, 
via (gHJ), to 

d ( i)3& = (84) 

and by taking the trace we obtain 

<5 (2) sft = -d (1) 3ft c (85) 
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(or equivalently using the symmetry of the Riemann tensor in ([83)1 ). 

By following a procedure similar to that leading to ()80|) exploiting the de Rham oper- 
ator for tensor valued 2- forms [HI HE IIH H] , which in our notation is simply A( 2 ) (or A^)), 
by using (|82l85j) . one easily derives 

A (2 )»= -d (2 )d(i)»ic. (86) 

For those spaces with V[ a Rb] c = 0, 

d ( i)3fhc = A ( 2)3? = 0. (87) 

This Laplace-like equation for the Riemann tensor appears to have been written down 
first of all by Penrose [17], motivated by spinors, and later in jUj; in both cases, only four 
dimensional Ricci-flat spaces were considered. Subsequently the general version has 
been presented in jSJ HHl EH1 El , as a direct result of differentiating the Bianchi equations. 
The typical index version of (|86j) reads 

V e V e R ab cd + 4R^ ef[c R b]ef d] + R abef Re fcd - 2R ab e[c R e d] = 4R la [c b] d ] • (88) 

(One should note different sign conventions, and some possible sign inconsistencies in the 
literature; see pQ for a summary.) 

As explained earlier in Sectional the trace and index symmetry properties across the 
two pairs of indices do not commute with the operator A/ 2 ); and so in particular, 

(z) (A (2) ft) oM + 0, {%%) tr(A (2) &) + A (2) tr(&) = A (2) 3foc 

as is easily confirmed. However, this operator A( 2 ) has been proposed in [HJ 01 [HI [5| as 
the natural analogue for the Riemann tensor of the usual de Rham operator A for the 
electromagnetic field tensor; we would argue, precisely because of these deficiencies, that 
it is not the most appropriate operator. It should be understood that in these references 
only vacuum spaces are being explicitly considered, and in such spaces the commutativity 
properties (ii) which fail in general are satisfied trivially because of the absence of Ricci 
tensor terms. 

Of course, we could have done the whole calculation using the corresponding trans- 
poses, and we would easily arrive at 

A (1) K = -d {l) d {2) ^ic (89) 

which is strictly equivalent to ([86)1 . due to the symmetry of the Riemann tensor pairs. 

A similar argument could be given concerning the information obtainable using the 
weighted operator A. Despite this fact, in general we believe that the operators Am and 
A(2) are potentially confusing and not really particularly suitable; this is especially true 
for general tensors with no symmetry between blocks, but also for the Riemann tensor 
itself. We firmly believe that the proper operator to be used here is A. From (|%ol l%9"|) it 
follows immediately that 

Aft = - - (d(i)d( 2 ) + d(2)d(i))9fo"c . (90) 
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This is a more useful version of a Laplace-like equation for the Riemann tensor, and it 
does not have the disadvantages of (jHUj) . It is easy to check, for instance, that the two 
examples considered above lead immediately to trivial identities: 

(z) (Aft) 0[M = 0, [ii) tr(AK) = Atr(9?) = M.ic. 

It should be remarked that the second equation here implies, in fact, that there is no 
non-trivial expression for Ajfoc derivable from (|9U|) . If one wishes to find a Laplace-like 
equation for the Ricci tensor, the best route is to use (|37l 155)1 to take the trace of ()83|) to 
get 3 

5 {1) mc = ~d i2) R (91) 

where R = tr(5J) = -R c c is the scalar curvature; from which, together with it follows 
that 

A (1) ®zc = -5(i)5 (2) 5i - ^d {1) d {2) R. (92) 
For the same reasons as above we prefer to combine it with its A( 2 ) counterpart 

A( 2 )3foc = — 5(2)5(i)3fa'c - -d(2)d(!)R (93) 

to obtain 

A3fac = --(5(1)5(2) + 5(2)5 ( i))K - -(d(i)d (2 ) + d^d^R . (94) 

A final contraction here can be obtained by remembering that A commutes with contrac- 
tion, and using (|38jl on the first term on the right hand side; on the otherhand, we cannot 
use (pT7j) on the second term since R is a scalar — acting here as a double (0, 0)-form. 
However, in index notation [(d^d^ 2 ) + ^(2)^(1) )R]ab = R,ab + R,ba and hence its trace is 
2AR. Hence we obtain the usual Laplace-like equation for the scalar curvature R 

^AR = --(5(i)5(2) + 5 (2 )5(i))9foc = -5(i)5( 2 )3fa*c , (95) 

where the last equality is due to the symmetry of the Ricci tensor. 
The index version of (j9~U|) was written in JHJ 5 PP and reads 

V e V e R ab cd + 4rt a ef[c R b ^ d] + R abef R efcd + R ab e[c R e d] + R cd e[a R e b] (96) 

= 2rt a [c ^ d] + 2R [c ^, d ^ 



3 This formula (|91|l is an illustrative example of the possibility mentioned on a number of occasions 
in this paper that some multiple forms must be considered as having extra blocks with no indices, when 
combined with other such forms. In this case, despite R being a scalar, one has to take its "rf(2)" -derivative. 
This means that R in this equation acts as a 'double (0,0)-form\ The general rule for these situations is 
to define the form-structure number of an equation as the maximum one for the different terms involved, 
and then every term in such equation has to be considered accordingly. See also Eas. (|92l95|l . I|71|) for 
other examples. 
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while the index versions of (}9"4"j) and (|95|1 are the familiar equations: 




— Raefb^ = V e V fR £ (a ft) , 



We believe that the equation (J9U|) is the natural Riemann tensor analogue of (|81|). and of 
course we also have 



Let us comment on the implications of the different definitions of 'harmonic tensors' for 
the Riemann tensor. From (J87)) . (|97|) we see that the Riemann tensor is fully harmonic 
(ergo harmonic) in spaces with d^^ftic = ( < ^=^ > c/( 2 )9fac = 0), including manifolds with 
a parallel Ricci tensor and its particular cases of Ricci-flat and Einstein spaces. As noted 
above, it has sometimes been the case that the operator A(2) has been considered as the 
extension to arbitrary tensors of the de Rham operator for single forms; and hence for 
the special case of the Riemann tensor the condition A( 2 )9ft = has been identified as 
its harmonic condition [3J |U] . However, we want to stress the following two important 
remarks, which highlight the properties of A and support our claim that this is the right 
operator, giving the proper harmonic condition: 

• The class of pseudo-Riemannian manifolds with a harmonic Riemann tensor is larger 
than that containing a fully harmonic Riemann tensor. This is clear from (|86|). (|89|) 
and (j9"Uj) . as we only need the condition (d^d( 2 ) + d( 2 )d(i))3?zc = for the former 
case. (Of course, for compact without boundary proper Riemannian manifolds the 
fully harmonic and harmonic cases are equivalent, as follows from Theorem 14.31 ) 

• We also note that the most common Laplacian used for arbitrary tensors, and 
thought to be mathematically appropriate, has been the Lichnerowicz Laplacian A^. 
We have pointed out, in general, that tensors harmonic in the sense of A^T = 0, 
do not coincide with (fully) harmonic tensors; in particular, as follows from (J3U|) . 
the Riemann tensor of pseudo-Riemannian manifolds (compact or not), are not 
harmonic by using Al, not even for Ricci-flat cases. Probably this was the reason 
that some authors chose equation (JHUj) (or equivalently (jHHl)) to define the Laplace- 
like or wave equations for the Riemann tensor, letting aside the probably 'more 
natural' operator A L . 

We believe that these two remarks taken together with other advantages mentioned be- 
fore demonstrate that the weighted operator A is better suited than either of Au^, A L , 
and probably the right choice for defining harmonic tensors and studying the derived 
implications. 

Finally, we want to remark about yet another good property of A (shared in this case 
by A(!) and A( 2 )): for any double (2,2)-form 1Z satisfying the first Bianchi identity (that 




(97) 
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is, a Riemann candidate), its 'Weyl part' is uniquely defined as usual by taking out the 
traces. We can denote by jy{7£} the Weyl part of a Riemann candidate. For example, 
VTjK} = C where C denotes the Weyl tensor. As proved first in j^j, the operator W 
commutes with A( 2 ), and therefore with A. Hence, this provides the straightest route to 
the correct Laplace-like equation for the Weyl tensor, because 

W{AU} = A(W{^}) = AC 

so that, on using (jUU|) we obtain 

AC = --W {{d m d m + d {2) d {1) )mc} . 

Now, it is very easy to produce an index version of this equation. (See pQ for different 
versions) . 



8 Summary and Discussion 

The inspiration for this paper was the result — that there exists a single potential for 
any Weyl candidate in any dimension — which we obtained in [21] in a rather pragmatic 
manner; here, that result has been shown to be a special case of a much more general result, 
which itself is a consequence of a significant generalisation and innovation in formalism. 
The underlying approach has been the systematic consideration of tensors as r-fold forms 
which has been explained and discussed at length in [15] , 

The generalisation of the differential form approach to tensor-valued forms and the 
extension of the use of the exterior differential d to such quantities is well known; we 
have emphasised the extension of the use of 5 and A also, and highlighted how an re- 
fold form can be thought of as r different tensor-valued forms, by taking each of the 
r sets of antisymmetric indices as the form indices, and defining the three operators 
d(i), S(i) and A(j), i = l,...r, associated with each in turn. However, we believe the 
explicit introduction of the weighted de Rham operator A as defined in this paper adds 
an important new ingredient. Two crucial properties of this operator are that it enables 
us, via a simple Laplace-like equation, to define an associated superpotential of exactly 
the same tensor type, and in addition to define potentials in a very natural manner, 
other generalised Laplacian operators for tensor-valued forms lack one or both of these 
properties. Hence we believe that A is a more powerful operator than A^, and so will be 
very useful in the type of formal investigations where Al has been used previously [4*T] . 
|42| . jjj, and as a powerful alternative in investigations of harmonic tensors. 

As a matter of fact, we have shown that the harmonicity property derived from the 
use of Al is different from that derived from A, and that the latter has better proper- 
ties concerning the harmonicity of Riemann curvature tensors (Section EJ) and in proper 
Riemannian manifolds ( Theorem 14. 3|) . We have also proved that one can use the defin- 
ing properties of A to define potentials and thereby to construct a generalised Hodge 
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decomposition. Theorem 14.21 provides this fundamental generalisation to arbitrary ten- 
sors (viewed as r-fold forms) of the local Hodge decomposition for antisymmetric tensors 
(single p-forms); an obvious next question to be considered is the possibility of a gener- 
alisation of the global Hodge decomposition for closed proper Riemannian manifolds. A 
first partial step towards this direction has been established in Theorem 14.31 which we 
believe was not known hitherto. 

The remaining task will require, on the one hand, a systematic analysis of the powers 
of the operators du) and S^, and on the other hand, a complete resolution of the relation 
between tensors with the property d^T = 0, and those with the property 3S : T = d^S 
(and the dual analogue with instead of <%).) By the way, this will allow us to see 
if there are some implications of the Bianchi identities ()82I84|) for the potentials defined 
in Theorem 16.31 in the same way as the first of the Maxwell equations (|8U|) implies the 
existence of a unique electromagnetic potential such that F = dA. All these matters are 
closely related, actually, to the gauge problem for the general potentials that we have 
mentioned in the paper and is the subject of current investigation [2E]- We are going to 
comment along these lines now. 

As is well known, the operators d and 5 are clearly nilpotent for single p-forms in 
curved spaces, and for tensor-valued forms in flat spaces. Deeper nilpotent properties for 
tensor-valued forms in curved space, and deeper commutator relationships with A, and 
other operators (such as the first order operator which defines the potential for the Weyl 
tensor) also need to be explored. 

As an application of the general result for all tensors in Theorem 14 .2\ we considered the 
Riemann curvature tensor showing that it can be written in terms of a pair of potentials, 
as given in Theorem 16.31 (There had earlier been some suggestions that the Riemann 
tensor could be written in terms of one potential via the Riemann-Lanczos equations 10 , 
but subsequently this was shown not to be possible [201, " ~ a ^ l eas t f° r dimensions 
n > 3; the fact that such a single potential cannot exist is now better understood from this 
new result.) The usefulness of this pair of potentials, and in particular, the implications 
which the Bianchi equations impose on their relationship, need further study. 

Recently there have been attempts in three dimensions to obtain (2, l)-form potentials 
L a \j C for the Ricci tensor via the Ricci-Lanczos equations |13j . 



The reasons for the lack of success in obtaining such potentials for arbitrary Ricci tensors 
in become apparent when ()98|) is compared with ([71)1: from the latter it is clear that 
this type of potential formulation of the Ricci tensor requires two potentials. On the 
other hand, we know from (J72j) that the traceless part of the Ricci tensor can always 
be given in terms of one potential, in arbitrary dimensions; hence the partial success in 
[TB] in investigating special spaces in three dimensions with zero Ricci scalar. Moreover, 
we can see directly that, by excluding the Ricci scalar from the potential structure, we 
will avoid the global obstructions involving the Ricci scalar in |T3] i n all dimensions (and 
the local obstruction in |22] in three dimensions). It is also instructive to compare the 
complicated nature of the second order equations for the various superpotentials in [T3] 
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with the corresponding simple and natural counterpart (|40p involving the weighted de 
Rham operator A. 

There are many interesting aspects of the new potential for the Weyl tensor which 
require deeper analysis. When the new potential for the Weyl tensor in all dimensions 
was obtained in [27] we pointed out that, in four dimensions, this new (2, 3)-form potential 
P ab c de is identical to the double dual of the classical Lanczos (2, l)-form potential H ab c 
[4"Uj . [55] , [2I|, [21- It is now clear why attempts to find higher dimensional analogues 
of the Lanczos potential gave negative conclusions [2E], [23] : we should have been 
looking for higher dimensional analogues not of the Lanczos potential, but of its double 
dual. There do not seem to be deeper implications for this new result in four dimensions; 
investigations of the Lanczos potential in four dimensions seem to be most efficiently 
dealt with in spinors |38j . [2]. However there are very important implications in higher 
dimensions. 

As emphasised in j2TJ, we now have, in all dimensions, an explicit potential for the 
Weyl tensor which supplies a tensor which is an 'integral' of the Weyl tensor at the level of 
the connection, and whose 'square' has units L~ 2 which are precisely the units we would 
expect for gravitational 'energies'. This suggests the usefulness of the wave equation and 
the super-energy tensor of this potential, in all dimensions with Lorentz signature, as an 
alternative to the Bel-Robinson tensor for such mathematical investigations as positivity 
properties, stability and the Cauchy problem for the Einstein equations [121 El- 

We have built our results on second order Laplace-like equations, with appeals to 
the Cauchy-Kovalewski theorem J5] so that our results are local for analytic pseudo- 
Riemannian metrics. However, we expect that these results can be generalised in the 
usual way; from the point of view of general relativity we can appeal to stronger theorems 
[3Tj . [33] when we specialise to spaces with Lorentz signature. However, of course the 
potentials are by definition first order, and now having established their existence it 
would be more natural to consider them in a first order system. In particular, the single 
potential for the Weyl tensor is an attractive candidate for deeper analysis in this context; 
preliminary investigations indicate that more direct and powerful results can be obtained 
by treating this definition as part of a first order symmetric hyperbolic system. 

When we choose a potential we know that it is not unique, and a full understanding 
of the role of gauge, which is more complicated for tensor- valued forms than for single 
p-forms, will be very important for further work. In particular we will need to have a 
set of explicit gauge equations to complete the first order symmetric hyperbolic system 
for the Weyl tensor. Furthermore, for the Weyl tensor, we would hope for a second-order 
linear equation for its potential with principal part of type V V^P cdei so that this will 
give an elliptic equation for positive-definite metrics and a wave equation for Lorentzian 
signature; in order to obtain such a simple version of the second-order linear equation for 
the potential we will need to exploit the gauge freedom. In curved four dimensional space 
(where we can be guided by the Lanczos potential [SHI, [21 > [21]), and in n-dimensional flat 
space it is easy to identify this gauge freedom and obtain the expected simple Laplace- 
like equation for the potential, but for other dimensions in curved space the calculations 
become long and complicated [28] . 
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